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Abstract 



In this paper we develop mathematical models for collective cell 

motility. Initially we develop a model using a linear diffusion-advection 

O type equation and fit the parameters to data from cell motility assays. 

^ This approach is helpful in classifying the results of cell motility as- 

qh say experiments. In particular, this model can determine degrees of 

""^ directed versus undirected collective cell motility. Next we develop a 

^-H model using a nonlinear diffusion term that is able capture in a unified 

^ way directed and undirected collective cell motility. Finally we apply 

the nonlinear diffusion approach to a problem in tumor cell invasion, 
Q^ noting that neither chemotaxis or haptotaxis are present in the system 

"^ under consideration in this article. 

o 

S^ 1 Introduction 



A basic characteristic of many organisms is their capability for active move- 
ment known as motility. Cell motility is fundamental in many important 
C^ and interesting biological phenomena such as morphogenesis, wound healing, 

and tumor invasion and metastasis (3]-[5]. The association of cell motility 
with diseases, and in particular cancer, makes it of interest not just for basic 
science but also in medical research. 

Much of the experimental research into cell motility focuses on molecu- 
lar or mechanical details and is carried out in such a way that it is difficult 
to form a unified picture of cell motion |3|. To complicate matters further, 
both individual cell motion and collective cell motility are important in many 
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phenomena. For example, the conventional view of tumor metastasis is that 
individual cells detach from a primary tumor and migrate to other areas |4|. 
However, now it is known that sheets or clusters of cells that maintain cell-cell 
contacts and migrate collectively also influence tumor invasion and metas- 
tasis |4l[5]. Spurred by experimental work and the surrounding difliculties, 
researchers have begun to develop mathematical and computational models 
of cell motility (see (ll|3}[6] and references therein). Such models can aid 
experimentalists in hypothesis testing, experiment design, and the creation 
of a unifled picture of cell motility. A survey of mathematical and computa- 
tional models of cell motility can be found in the article p!o|. It would also 
be valuable if such models could be used to provide insight into systems of 
interest in biomedical research such as tumor invasion. 

Mathematical models of individual cell motion such as in [6, 10 are able 



to incorporate much of the mechanical details of migration and model the 
migration of cells in tissue in which motion in some directions may be re- 
stricted. For example, in [g] Hillen derives mathematical models for three 
dimensional cell motion through temporally varying tissue networks using 
kinetic transport equations. In this case, the motion in certain directions 
is limited by the orientation of tissue or extra cellular matrix. This type of 
motion is often referred to as mesenchymal motion, and direction of motility 
is influenced by external factors not the intrinsic behavior of the cell popu- 
lation. Two cases are studied in |6|, cell motion at the individual cell level 
and at the population level. Hillen begins with a transport equation model 
for individual cell mesenchymal motion and then uses scaling arguments to 
derive a model for motion of a population of cells. 

In this work we consider the development of models of collective cell 
motility that can be applied to the case of a group of cells that undergo 
collective migration where the degree of directed motility is determined by 
the intrinsic behavior of the cell population. More precisely, we consider the 
epithelial-mesenchymal transition (EMT) which is characterized by the loss 
of cell-cell adhesion and the resulting effects on motility [4|[7[[l7ll. The work 
presented in this paper is motivated in part by experiments of the Stipp 
lab [8[[9l pT|p!4] and in part by the possibility of applying the results to simple 
but useful models of tumor invasion such as in [T3[[T5|[T6]. We note that 
the cell motility scratch assays by the Stipp lab include neither chemotaxis 
since no chemoattractant is introduced into the scratch. Likewise, for the 



models in 13, 15, 16 of spherical in vitro tumors no chemotaxis or haptotaxis 



since no chemoattractants or macromolecules exist outside the boundary of 



the tumor. In particular, the model of Stein et al. [13] is developed for 
glioblastoma tumor spheroid in a three-dimensional in vitro experiment. 

2 Biological Background 

Cell adhesion, is the result of the interactions of certain cell receptors with ad- 
hesion molecules. For example, proteins know as cadherins have been shown 
to influence cell-cell adhesion of cancer cells (4J . In particular epithelial cells 
express the protein E-cadherin which is responsible for the tight junctions 
between cells that is characteristic of the epithelium. Collective cell migra- 
tion depends on cell-cell adhesion, while migration of individual cells is often 
associated with the loss of cell-cell adhesion. Such a loss of cell-cell adhesion 
may be an eflFect of an epithelial-mesenchymal transition (EMT). "An EMT is 
a culmination of protein modiflcation and transcriptional events in response 
to a deflned set of extracellular stimuli leading to a long term, albeit some- 
times reversible, cellular change. Core elements of EMT include reduction 
of cell-cell adherence via the transcriptional repression and delocalization of 
cadherins (adherens junctions), occludin and claudin (tight junctions), and 
desmoplakin (desmosomes)." |7|. It is thought that in (epithelial) cancers an 
EMT is essential for invasion and metastasis. 

The Stipp lab has performed experiments known as scratch assays in 
order to study the effects of an EMT on collective motility properties of 



cells 14 . In a scratch assay, cells are cultured on a surface and then a 
scratch is made to create a gap between some of the cells. Cells are then 
filmed migrating to "fill in the gap". Software (NIH Image J 1.63 software) 
can then be used to track some number of the cells in order to generate data 
from which cell velocities can be determined [T4|. That is, the experiments 
provide data about the individual random walks of the cells. By applying 
standard statistical techniques to this data, one can compute diffusion and 
drift coefficients that offer a potential method for characterizing the results 
of such scratch assays. In particular, by computing the drift and diffusion 
coeflBcients we may be able to characterize the degree of cell-cell contact in a 
given assay. This is based on the observation that cell-cell contact infiuences 
motility properties (i.e. how quickly a population of cells migrates to fill 
in the gap). Using the mathematical models described below the motility 
properties can be characterized. 

Figure [T] shows frames from a film of a scratch assay performed by the 




Figure 1: Scratch assay. Figure shows stills at labeled times from a film 
of a scratch assay. Stills in the left hand column show PC-3 (epithelial) cells 
while the right column shows TEM4-18 (mesenchymal) cells. Note that cells 
on the right show reduced cell-cell adhesion. Figure courtesy of Christopher 
Stipp. 



Stipp lab. "The cell types were PC-3 (epithehal) and TEM4-18 (mesenchy- 
mal). TEM4-18 cells are a subline of PC-3 cells that was obtained experi- 
mentally by collecting a subpopulation of PC-3 cells with enhanced ability 
to traverse an endothelial cell layer. Later it was found that TEM4-18 cells 
had a gene expression profile consistent with a loss of epithelial cell charac- 
teristics, such as stronger cell-cell contacts, as compared with the parental 
PC-3 cell population as a whole. tl4|" The creation of the TEM4-18 subline 
from PC-3 cells is discussed in 191. 



3 Mathematical Modeling 

Using the data from the scratch assays, we have developed mathematical 
representations of collective motility. One such representation simply makes 
use of the diffusion and drift coefficients by input into a linear equation of 
the form 



dtU = V • (DVu-vu), 



(1) 



where u{x^t) represents the cell density at a given position x and time t, 
D is the diffusion coefl&cient, and v is the drift velocity. An alternative 
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Figure 2: Results of the drift-diffusion model of cell motility. This 
figure shows the front propagation of invading cells based on model ^ and 
data from scratch assays. The front propagation characterizes the motility 
properties of cells in a scratch assay Since motility properties are linked with 
cell-cell contact, this allows one to characterize the degree of cell-cell contact 
by the drift and diffusion coefficient determined from the assay 



approach developed in [12j uses a discrete modeling method to build cell 
motility models based on individual cell random walk data. 

Figure [2] shows the results typical of using ^ to track the front propaga- 
tion in a scratch assay. Here we have used a scaling to have unit maximum 
cell population density. The spatial domain for the computation is the in- 
terval [—15,1]. On the left boundary we take u{—15^t) = 1 for all t, a 
so-called far-field condition. This is justified, since in a typical scratch assay 
the scratch is small compared with the area containing cells. We employ a 
no-fiux condition on the right boundary. We take as the initial condition 



u{x^ 0) 



ifx<|. 



1, 

0, otherwise. 



(2) 



This represents the assay at the time of the scratch. Here we only consider 
cells to the left of a scratch as in figure [TJ This is sufficient due to the left- 
right symmetry about the center of the scratch. The scratch creates a front 
between a region with cells and a region without cells, which moves as the 
cells migrate to fill in the gap. Figure [2] shows a right-moving cell population 
that corresponds to the motion of the front in the scratch assays. We obtain 
from the computations the front speed and compare it with experimental 



results. In this case parameter values are chosen to correspond to a scratch 
assay with epithelial (PC-3) cells. This provides a method for characterizing 
the motility properties, and possibly the degree of cell-cell adhesion, in such 
assays. 

In [2] , Aronson describes the role of diffusion equations in modeling dis- 
persing populations. In particular, [2] covers how the behavior of a dispersing 
population is related to the form of the diffusion operator appearing in the 
equations modeling the spatial dynamics of a population. It is pointed out 
that it can be useful to consider nonlinear diffusion terms. In |1| the authors 
present a discrete model of cell motility that incorporates cell-cell adhesion. 
This discrete model has as its continuum limit a nonlinear diffusion equation 
of the form 

dtU = V •{D{u)Vu). (3) 

Due to the occurrence of cell-cell adhesion, in |1|, the diffusivity is described 
by a function of the form 

D{u) = a{u-bf + c, (4) 

where a, 6, c are constants related to the adhesion properties of cells. Under 
the influence of fll[2] we seek an alternative to (fTl) by considering nonlinear 
diffusion equations of the form ^ . 

Figure [3] shows the results of using a nonlinear diffusion model of the 
form ([3]) to represent cell motility as observed in the scratch assays discussed 
above. In order to get a qualitative match with the results, such as in figure 
[21 using (fTl) we have taken the nonlinear diffusivity to have the form 

D{u) = a • max(0, u — b)u + c, (5) 

which is similar to ([4]). The boundary and initial conditions are the same 
as above. Figure |3] again simulates a right-moving cell population (epithelial 
(PC-3) cells) that corresponds to the motion of the front in the scratch assays. 
The advantage of using a nonlinear diffusion model is that the effects of the 
degree of cell-cell adhesion on motility, of interest in the study of the EMT, 
can be explicitly represented. Note that by making the constant a smaller 
or the constant b larger in ^ the diffusion becomes closer to linear, that is 
we get closer to purely random motility. 




Figure 3: Results of nonlinear model of cell motility. This figure shows 
the front propagation of invading ceUs based on model ^. 

4 Application to Tumor Invasion 

In [l5|[T6J Swanson et al. apply simple but useful mathematical models of 
tumor invasion. In particular the authors use a reaction-diflFusion equation 
of the form 



dtU = V ' (DVu) +gu[l 



u 



Ur, 



(6) 



to model the invasion of a spherically symmetric gliobastoma tumor. The 
use of a reaction-diffusion equation is based on the assumption that the 
net dispersal of tumor cells is captured by a random walk which translates 



mathematically into diflFusion. In the case of [15, 16 diflFusion is assumed to 
be linear. The term gu{l ^^—) describes the proliferation of cells. Based 



on invasion assays of tumor cells, Stein et. al. 13 suggest that there is a 
need to consider not only random motility of invasive cells but also directed 
motility. In particular the authors of |13j conclude that while diffusion alone 
is sufficient for describing cell motility of cells making up the central core 
of the tumor, cells at the invasive rim of the tumor have different motility 
properties that include both random and directed motility. By modifying ^ 
the authors of fls] arrive at the following model for the invasion of tumor 
cells 



dfU = V • {DVu — vu) + s5{x — xq — vt) + gu i 1 



u 



u^ 



(7) 



Here 6{-) is the Dirac delta function. In the above model the s6{x — Xq — vt) 
term represents the shedding of invasive cells from the invasive rim. 

Figure [4] shows simulations using ([7|). Here we have assumed a spherically 



symmetric tumor as done in [13, 15, 16). The interesting feature of these 



results is the behavior of the invasive rim. We note that in |13| the cells 
forming the tumor core are treated separately, as a result the behavior of the 
cells making up the core are not shown in figure 111 In |13l , the parameters 
D,v, 5,5' in ([7| are optimized to fit experimental results obtained by the 
authors. 

Based on the scratch assay experiments and the nonlinear model for cell 
motility presented above we modify the model of Stein et al. ([7| to have the 
form 

dtu = \/ ' {D{u)\/u) + s5{x -xo-\t)+gull — j , (8) 

where the motility is described by a nonlinear diffusion term of the general 
form V • {D{u)Vu). We interpret this as having the motility dependent 
on the population density. This is of particular interest in the case where 
there is cell-to-cell adhesion. Figure [5] shows simulations using the model 
([s]). We note the qualitative similarity of the behavior of the invasive rim. 
One difference however, is that we do not treat the cells at the tumor core 
separately. This is the reason for the difference in centers of the simulations. 
Figure [5] does not have the "hole" in the center that figure [4] does. 

In conclusion, based on the data from the scratch assays performed in the 
Stipp Lab we are able to obtain information on the random walks of individ- 
ual cells. From the individual random walks we are able to classify the results 
of the experiments in terms of the diffusion and drift coefficients. Using these 
coefficients we can represent the motility of cells by a linear drift-diffusion 
model. This model is capable of providing a qualitative match to experimen- 
tal observations. Further, by using a nonlinear diffusion model, we can again 
represent the motility of cells such that we again obtain a qualitative match 
to experimental observations. The value of the nonlinear model is that it 
includes motility that is population density dependent, which is relevant in 
the case that there is cell-cell adhesion. Moreover, for certain applications 
such as tumor invasion the nonlinear model may be appropriate since as 
discussed in |5] cell-cell adhesion and hence population density dependent 
migration occurs. The models developed here provide an alternative for cap- 
turing directed cell motility in cases where mechanisms such as chemotaxis, 
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Figure 4: Simulation of tumor invasion with Stein et al. model. This 
figure shows simulations (using ([7|) of invasion of a sphericaUy symmetric 
tumor. The simulations are meant to model experimental results described 



in 13 . The parameter values from |13| are used. 
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Figure 5: Simulation of tumor invasion with nonlinear model. This 
figure shows simulations (using ([s])) of invasion of a sphericaUy symmetric 
tumor. The parameter values are chosen to match qualitatively with the 
behavior of the invasive rim in figure [4] 
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haptotaxis, etc. are not dominant, and when there are no external factors 
such as anisotropic tissue to influence the direction of ceU motion. 
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